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Using the 0{a) Symanzik improved action an estimate is given for the strange quark mass for unquenched 
(n/ = 2) QCD. The determination is via the axial Ward identity (AWI) and includes a non-perturbative evaluation 
of the renormalisation constant. Numerical results have been obtained at several lattice spacings, enabling the 
continuum limit to be taken. Our results indicate a value for the strange quark mass (in the MS'-scheme at a 
scale of 2 GeV) in the range 100 - 130 MeV. A comparison is also made with other recent lattice determinations 
of the strange quark mass using dynamical sea quarks. 



1. INTRODUCTION 

Lattice methods allow, in principle, the com- 
plete 'ab initio' calculation of the fundamental pa- 
rameters of QCD, such as quark masses. However 
quarks are not directly observable, being confined 
in hadrons and are thus not asymptotic states. So 
to determine their mass necessitates the use of a 
non-perturbative approach - such as lattice QCD 
or QCD sum rules. In this brief article, we report 
on our recent results for the strange quark mass 
for 2-flavour QCD in the MS'-scheme at a scale of 
2 GeV, mf^(2 GeV) (further details can be found 
in CI). 

The present phenomenological status is sum- 
marised by the Particle Data Group in [2| giv- 
ing an estimate for the strange quark mass of 
80 MeV < mp(2GeV) < 130 MeV. This is a 
large band, and it is be hoped that lattice compu- 
tations will reduce this significantly in the coming 
years. 



2. RENORMALISATION GROUP 
INVARIANTS 

Being confined, the mass of the quark, m^(M), 
needs to be defined by giving a scheme, <S and 
scale M, 



BABE 
q ■ 



(1) 



and thus we need to find both the bare quark 
mass and the renormalisation constant. An added 
complication is that the MS'-scheme is a pertur- 
bative scheme, while more natural schemes which 
allow a non-perturbative definition of the renor- 
malisation constants have to be used. It is thus 
convenient to first define a (non-unique) renor- 
malisation group invariant (RGI) object, which 
is both scale and scheme independent by 



where we have 



(2) 



*Talk given by R. Horsley at the Workshop on Computa- 
tional Hadron Physics, (Nicosia, Cyprus, September '05). 
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Figure 1. One-, two-, three- and four-loop results for 
[AZ*p(/i)]-i in units of A*^. 

The f3^ and functions (with leading coeffi- 
cients —bo, dmo respectively) are known perturba- 
tively up to a certain order^. In the MS scheme 
the first four coefficients are known, [311, and 
this is also true for the RI'-MOM scheme |5I()| 
(which is a suitable scheme for both perturbative 
and non-perturbative, NP, applications, see sec- 
tion E)). Note that in the Rl'-MOM scheme we 
also choose to expand the (3^ and 7^ functions 
in terms of g"^, [HI; other choices, of course, are 
also possible. 

In Figs. HandUl we show the results of solv- 
ing eq. as a function of the scale M = and 
M = fip for both the JTS and Rl'-MOM schemes 
respectively. We hope to use these (perturbative) 
results in a region where perturbation theory has 
converged. 2GeV corresponds to /Lt/A"^ ~ 8, 
where it would appear that the expansion for 
the M^-scheme has converged; for the Rl'-MOM 
scheme using a higher scale is safer (which is cho- 
sen in practice). However, when the RGI quantity 
has been determined we can then easily change 
from one scheme to another. Of course these 

^Analgous definitions also hold for other operators, see 
section |3 



Figure 2. One-, two-, three- and four-loop results for 
[AZ™'-™"(/xp)]-i in units of A*^. 

scales are in units of A*^^ which is awkward to 
use: the standard 'unit' nowadays is the force 
scale tq. To convert to this unit, we use the re- 
sult for tqA"^ as given in [7j. Of course, we must 
also give the physical scale. A popular choice is 
To = 0.5 fm, but there is some variation in possi- 
ble values, see section |S1 

3. CHIRAL PERTURBATION THEORY 

We have generated results for n/ = 2 degener- 
ate sea quarks, together with a range of valence 
quark masses. Chiral perturbation theory, xPT, 
has been developed for this case, |8I9| . We have 
manipulated the structural form of this equation 
to give an ansatz of the form 

romf = 

c^' [(^'0^^+)^ + {romiiof - (rom^+)^] + 

(rom^+)2 + 
i(cr + cD [irom^.f + (rom^o)2]' - 
{c^ + c-f'){rom^^r - 

[[rarriK+T + {roniKaf - (rom„+)^] x 
In {ironiK+Y + (^0™Ko)^ - {rom^+Y) + 
cf'{rom^+f\n{rom^+f, (4) 
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where 



en 

cf {{raml^f - 2{rampsf) ln(romp,)2 . (5) 



TUps, rripg are the valence and sea pseudoscalar 
masses respectively (both using mass degenerate 
quarks). The first term is the leading order, LO, 
result in xPT while the remaining terms come 
from the next non-leading order, NLO, in xPT. 
Thus we see from eq. (O that to NLO, we can first 
determine c^' and cf^' , i = b,c,d using pseu- 
doscalar mass degenerate quarks and then simply 
substitute them in eq. 

4. THE LATTICE APPROACH 

Approaches to determining the quark mass on 
the lattice are to use the vector Ward identity, 
VWI (see e.g. JOl)' where the bare quark mass is 
given in terms of the hopping parameter by^ 



1 

2a 



(6) 



or the axial Ward identity, AWI, which is the ap- 
proach employed here. Imposing the AWI on the 
lattice for mass degenerate quarks, we have 



(7) 



and A and V are the 0{a) improved^ unrenor- 
malised axial current and pseudoscalar density 
respectively and fhq is the AWI quark mass. So 
by forming two-point correlation functions with 
V in the usual way, this bare quark mass can be 
determined 



anio 



2{v{t)vm 



(8) 
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^This is valid for both valence and sea quarks, k^^ 
defined for fixed f} by the vanishing of the pseudoscalar 
mass, i.e. mps{K^^, k^^) 



■ 0. Kgj, has been determined in 



^The improvement term to the axial current, d^P together 
with improvement coefficient ca, (llj has been included. 
The mass improvement terms, together with their associ- 
ated difference in improvement coefficients, 6^, bp appear 
to be small and have been ignored here. 



We have found results for four /3-values: 5.20, 
5.25, 5.29, 5.40, each with several (three or more) 
sea quark masses and a variety of valence quark 
masses, pQ. 

Furthermore upon renormalisation we have 

A^ = ZaA^ , V^iM) = Zf,{M)V , (9) 

giving from eqs. © and 

Za 



zr = ^zUm) 



Zf{M) 



(10) 



As mentioned before, we use the Rl'-MOM 
scheme, H]. This scheme considers amputated 
Green's functions (practically in the Landau 
gauge) with an appropriate operator insertion, 
here either A or P. The renormalisation point 
is fixed at some momentum scale = /Zy, and 
thus we have 
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where Tq are one-particle irreducible (IPI) vertex 
functions, and Zq is the wave-function renormal- 
isation. (Our variation of the implementation of 
this method is described in 52] •) This determines 
Za and Zp ^ from which a chiral extrapola- 
tion, here using the sea quarks only, may be made 
to the chiral limit. For Za we make a linear ex- 
trapolation in aruq, Za = Aa + BAairiq, while for 
^m'-MOM must first subtract out a pole in the 



quark mass, [T5| 
metry breaking. 

MOM \ — 1 



which occurs due to chiral sym- 
We thus make a fit of the form 



) — Ap + Bp/airiq. We now have all 



the components, namely A/ 



zm'-MOM 



necessary to compute Zf^' and hence rom^' 



Ap and AZ? 

In 

Fig. El we show Z^' for /3 = 5.20, 5.25, 5.29 and 
5.40 . These should be independent of the scale 
afip at least for larger values. This seems to be 
the case, we make a phenomenological fit to ac- 
count for residual effects. 

5. COMPARISON OF Z?" WITH 
OTHER METHODS 

As many computations of the strange quark 
mass have used tadpole improved perturbation 
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Figure 3. Z^' for /3 = 5.20 (filled circles), /3 = 5.25 (filled squares), /3 = 5.29 (filled upper triangles), f3 = 5.40 (filled 
lower triangles) together with fits F(afj,p) = ri + r2(a/ip)^ + r:j/{a^p)'^ . 



theory together with a boosted couphng constant 
for the determination of the renormahsation con- 
stant, it is of interest to compare our results ob- 
tained in the previous section with this approach. 
Our variation of this method, tadpole-improved 
renormalisation-group-improved boosted pertur- 
bation theory or TRB-PT, is described in |14|. 
Regarding the lattice as a 'scheme', then from 
eq. we can write 



AZr(a)m,(a), 



(12) 



where the renormalisation-group-improved 
AZf^{a) is given by eq. |(2Jl. Furthermore in this 
'lattice' scheme, we choose to use ~ 5o/"oc 
where ^ (|TrC/°} ([/° being the product 
of links around an elementary plaquette) rather 
than go, as series expansions in ga are believed 
to have better convergence. This is boosted 
perturbation theory. (We shall use chirally ex- 
trapolated plaquette values as determined in ^ 
at our /3-values and so we add a subscript 'c' 
to Uq.) In the tadpole- improved method, not- 
ing that renormalisation constants for operators 
with no derivatives are ~ uqc, which indicates 



that Z^'uq^ will converge faster then Z^' alone 
we re-write eq. © in the two loop approximation 



^egi{trb-ft) 



(13) 



1 , ^1 2 

bo 



1 91 



where qi = {bod^J^ - 6id™o)/(2&o&i) + 
(pi/4) (60/61) with pi = i being the first coeffi- 
cient in the expansion of uqc- d^ff may be found 
by relating the (known) perturbative result for 
Z^ to AZ^-. 

In Fig. 0] we plot Z^' versus (3. Our NP re- 
sults from section 0] are shown as filled circles. 
They are to be compared with the TRB-PT re- 
sults denoted by empty squares. While there is 
a difference between the results, it is decreasing 
and thus may be primarily due to remnant 0{a^) 
effects, which disappear in the continuum limit. 
That various determinations of Z^' have differ- 
ent numerical values can be seen from the results 



*The TRB-PT subscript in brackets is there only to dis- 
tinguish the results from those obtained in section Bl 
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Figure 4. 

^j?G/ versus 3. The black circles are the 
results from section|3 while the open squares are the TRB- 
PT results. Furthermore the open diamonds and triangles 
are the NP results from 1111 , using the two different results 
for the axial renormalisation constant, )15l . (The empty 
triangle results have been slightly displaced for clarity.) 

of (open diamonds and triangles). In these 
results two different definitions of the axial renor- 
malisation constant have been used, |15j. Za is 
computed when dropping certain disconnected di- 
agrams, while includes them. (The differ- 
ence between the two definitions is an 0[a}) ef- 
fect.) Using in Z^' leads, perhaps coinci- 
dently, to very similar results to our NP results. 

Investigating the possibility of O(a^) differ- 
ences a little further, we note that if we have two 
definitions of Z'^' then if both are equally valid, 
forming the ratio should yield 

yBGl{X) 

Rl^^^ = ^ + 0{a'). (14) 

where Z^' is the result of section0]and X is some 
alternative definition (i.e. TRB-PT, ALPHA- Z^, 
ALPHA- ). In Fig. we plot this ratio 
for these alternative definitions. The ro/a values 
used for the x-axis are found by extrapolating the 
ro/a results to the chiral limit. This extrapola- 
tion and results (for (ro/a)c) are given in p. 

We see that (roughly) all three ratios extrapo- 
late to 1 which implies that any of the four deter- 
minations of Z^' may be used. This includes the 
TRB-PT result. Of course other TI determina- 



Figure 5. versus (a/ro)l for X = TRB-PT (open 

squares), X = ALPHA- (open diamonds) and X = 
ALPHA- Z^O" (open triangles). 

tions might not have this property, and also their 
validity always has to be checked against a NP de- 
termination, so this result here is of limited use. 
It is also to be noted that different determinations 
can have rather different O(o^) corrections, so a 
continuum extrapolation is always necessary. 

6. RESULTS 

Armed with Z^' , we can now find rom^' and 
hence the ratio rom^' / (roirips)'^, using the values 
of ro/a given in [7j . In Fig. we plot this ratio 
(against (rorrips)^) for f3 — 5.29. Using eq. (0} to 
eliminate c^' in eq. (j^l in favour of c^' where 



(rom;f+)2 -I- (romxo)2 - {rom^+)^ ' 

gives rgTOj, directly^ to NLO in our fit function. 

We have restricted the quark masses to lie 
in the range (roWps)^ < 5, which translates to 
irips < 880 MeV, which is hopefully within the 
range of validity of low order xPT results. (Us- 
ing ro/a, rather than their chirally extrapolated 
values for example, tends to give less variation in 
the ratio rom^' /{rompsY so we expect LO xPT 
to be markedly dominant.) Varying this range 

^This is preferable to first determining and c'^' , i = 
b, c, d by using eq. J3 S'^id then substituting in eq. as 
the direct fit reduces the final error bar on rom^' . 
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Figure 6. rom^^ / (romps)^ against (rompa)^, together with a fit using eq. j^J for /3 = 5.29. Filled points represent 
valence quark results while unfilled points are the sea quark results. The dashed line (labelled '\/2mfc') represents a fictitious 
particle composed of two strange quarks, which at LO xPT given from eq. J3} by (r"o"^x+ )^ + i^O'm-K'^)'^ ~ (''0"^,r+)'^> 
while the dashed-dotted line (labelled 'm7r') representing a pion with mass degenerate u/d quark is given by rom^+ . 



from 4 to 6 and higher gave some idea of possible 
systematic errors. Thus finally, for each /3- value 
we have determined r^mf^' and can now perform 
the last extrapolation to the continuum limit. 

Our derivation so far, although needing a sec- 
ondary quantity such as r^/a for a unit, de- 
pends only on lattice quantities. Only at the last 
stage, with our direct fit did we need to give a 
physical scale to this unit. A popular choice is 
To = 0.5 fm. However there are some uncertain- 
ties in this value; our derivation using the nu- 
cleon gave rp = 0.467 fm and so to give some idea 
of scale uncertainties, we shall also consider this 
value. (The main change when changing the scale 
comes from the tqS in eq. 10} , as mf^' oc tq, while 
changes in AZ^^ are only logarithmic.) 

Using the value for [AZ,^(2 GeV)]-i obtained 
in section |21 from Fig. ^ to convert raf^' to 
m"^(2GeV) gives the results shown in Fig. [7| 
Also shown is an extrapolation to continuum 
limit. We finally obtain the result 

mf^(2 GeV) = (16) 



r 117(6) (4) (6) MeV for ro = 0.5fm 

\ lll(6)(4)(6)MeV for tq = 0.467 fm ' 

where the first error is statistical. The second er- 
ror is systematic ~ 3 MeV estimated by varying 
the fit interval for {rQWipsY ■ We take a further 
systematic error on these results as being covered 
by the different tq values of about ~ 6 MeV. This 
is to be compared to our previous result using 
the VWI, J^, which gave results of 126(5) MeV, 
119(5) MeV for tq = 0.5 fm and 0.467 fm respec- 
tively. These results and extrapolation are also 
shown in Fig. [7| 

7. COMPARISONS 

It is also useful to compare our results with 
the results from other groups. In Fig. [S] we show 
some results for nj = 2 and Uf — 2 + 1 flavours 
(keeping the aspect ratio approximately the same 
as in Fig. [7| A variety of actions, renormalisa- 
tions, units and scales have been used (so the re- 
sults have been plotted in physical units using 
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Figure 8. Results for m^^{2 GeV) versus fm^ using the AWI (upper plot) and VWI (lower plot) methods. The results 
are presented with the collaborations preferred units and scales. Circles (together with a linear continuum extrapolation) are 
from this work and IIUI : diamonds from lltil : squares from 1171 : up triangles from 1181 : down triangles from 1191 : left triangles 
from 1201 : right triangles from NPR denotes non-perturbative renormalisation, while TB-PT denotes tadpole-improved 
boosted perturbation theory. I20I21I are for rif =2 + 1 flavours; the other results arc all for nj = 2 flavours. 



the authors prefered values). In particular the 
HPQCD-MILC-UKQCD ^ and HPQCD ^ 
collaborations use improved staggered fermions. 
These fermions having a (remnant) chiral symme- 
try are in the same situation as overlap/domain 
wall fermions where there is no distinction be- 
tween VWI and AWI quark masses; the bare 
quark mass in the Lagrangian simply needs to 
be renormalised. 

As seen earlier in section |51 it is noticeable that 
the (tadpole improved) perturbative results lie 
lower than the non-perturbatively renormalised 
results. Also results with a < 0.09 fm (i.e. 

< 0.008 fm^) appear to be reasonably consis- 
tent with each other (this is more pronounced for 
the AWI results than for the VWI results). While 
results for a < 0.09 show some lattice discretisa- 
tion effects, using results at larger lattice spac- 
ings seems to give a fairly constant extrapolation 
to the continuum limit. A similar effect has also 
been seen elsewhere, for example in the determi- 
nation of tqA"^ for 71/ = flavours, 7 , where 



coarse lattices also show this characteristic flat- 
tening of the data. 

Finally, we compare these numbers with re- 
sults from the QCD sum rule approach. A re- 
cent review of results from this method is given 
in 122], citing as a final result m*^®(2GeV) = 
99(28) MeV. This covers the lattice results in 
Fig.EI 

8. CONCLUSIONS 

In this article we have estimated the strange 
quark mass for 2-flavour QCD and found the 
result in eq. (|16|l . using 0(a)-improved clover 
fermions and taking into consideration non- 
perturbative (NP) renormalisation, the contin- 
uum extrapolation of the lattice results and the 
use of chiral perturbation theory. The NLO chiral 
perturbation theory yields a correction of about 
5% to the LO result, and the relevant low energy 
constants are in rough agreement with the phe- 
nomenological values. 
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Fi gure 7. Results for mg^^(2GeV) (filled circles) versus 
the chirally extrapolated values of (a/ro)^ (as given in U]) 
together with a linear extrapolations to the continuum 
limit. For comparison, we also give our previous result 
using the VWI, IIUI (open squares). 

In conclusion, although there is a spread of re- 
sults, it would seem that the unquenched strange 
quark mass determined here is not lighter than 
the quenched strange quark mass and lies in the 
range of 100 - 130 MeV. 
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